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Abstract. - We consider a quantum quench in a system of free bosons, starting from a thermal 
initial state. As in the case where the system is initially in the ground state, any finite subsystem 
eventually reaches a stationary thermal state with a momentum-dependent effective temperature. 
We find that this can, in some cases, even be lower than the initial temperature. We also study 
lattice effects and discuss more general types of quenches. 
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• Introduction. — The problem of a quantum quench, 
i'.e. an instantaneous change in the parameters that 
determine the dynamics of a quantum system, has been 
recently studied in various theoretical models. These 
include systems of free bosons or fermions in which the 
cjuenched parameter is the energy gap [1-3], integrable 
models described by conformal field theory [1, 2, 4-6] 
and exactly solvable spin chains [4,7,8]. Several studies 
have also been carried out in non-integrable models using 
ilumerical methods or approximation schemes [9-12]. 
Perhaps one of the most interesting results is that in the 
thermodynamic limit it is possible under quite general con- 
ditions that for large times an arbitrarily large subsystem 
tends to a stationary state with thermal characteristics, 
a process we can call thermalization [1-3,13]. 

\ To be more specific, let us imagine a system of coupled 
harmonic oscillators or equivalently a free field theory, 
described by a general dispersion relation with some 
energy gap or "mass" mo and maximum group velocity of 
excitations c. Assume that the system initially lies on the 
ground state of the initial hamiltonian Hq and at time 
t — Q the mass is suddenly changed from mo to a different 
value m. For t > 0, the state of the system evolves 
according to quantum mechanics {unitary evolution), 
i.e. it is fully isolated from its environment. After this 
quench there is an extensive energy excess in comparison 
with the ground state of the final hamiltonian H, which 
is distributed over the excitation levels of H . It then 
turns out that the two-point correlation function, which 
in free systems contains all the information required 
to determine their state, acquires for sufficiently large 



times the form of the correlation function of a system 
at thermal equilibrium with a momentum- dependent 
effective temperature (3cs{k) [1, 2, 14]. Note that this 
momentum dependence is expected for a free theory since 
each momentum mode-excitation evolves independently 
from the others and can therefore thermalizc to a differ- 
ent temperature. For any local observable, the quantum 
interference between all the momentum modes gives rise 
to a single effective temperature, which, however, may be 
dependent on the chosen observable. However, the large 
distance behaviour of the two-point correlation function 
is determined only by the zero- momentum mode /3eff(0). 
In the so called deep quench limit where toq ^ 'ti, the 
effective temperature P~^{k) is asymptotically equal to 
mo/4 and becomes momentum independent [1,2]. 

In most studies made so far, the initial state of the 
system is the ground state of Hq. Here instead, we 
assume that, before the quench, the system is prepared in 
a thermal state. Then we perform the quench exactly as 
before keeping the system isolated from the environment 
(the couphng to the environment is studied in [15]). For 
brevity we can call this a thermal quantum quench, in 
contrast to the previous case that we will call a pure 
quantum quench. Few other studies including the effect 
of a thermal initial state exist [8, 12, 16], but consider 
different kinds of quenches. The solution of this problem 
is in general prohibitively difficult and effective numer- 
ical methods must be used to extract useful physical 
quantities. However, in free field theory the solution is 
particularly simple and allows us to investigate the effect 
of the initial temperature (3o on the time evolution of the 
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system. Especially we can study the crossover from the 
"/loi" limit, i.e. when the initial temperature dominates, 
to the "co^d" limit when it can be ignored, in which case 
we recover the pure quantum quench results. 

It turns out that after such a thermal quench the sys- 
tem also thermalizes in the sense described earlier. An 
intuitive guess would perhaps be that the final effective 
temperature /3~g is just the sum of the initial tempera- 
ture (3q^ and that corresponding to an equal pure quench 
of the mass However this is not true and l3~f^ can 

instead be lower than Pq^. In the hot deep quench limit 
{00^ ^ iTiQ ^ m), for example, the final temperature is 
shown to be half of the initial one. 

In the present paper we provide the solution of the prob- 
lem of a thermal quantum quench in a system of free 
bosons. We start with a simple harmonic oscillator where 
we calculate the analogue of the correlation function as a 
function of time. Then we use this to find the correlation 
function in a bosonic free field theory with a relativistic 
dispersion relation and study its large time as well as real 
space behaviour. We thus identify its stationary part and, 
by comparison to the thermal expression, derive the ef- 
fective temperature, whose properties in several limits are 
discussed. Furthermore we investigate several other types 
of quenches and derive general expressions for the mo- 
mentum distribution of excitations in lattice models also 
extending to the fermionic case. 

Harmonic oscillator. — We first consider a single 
harmonic oscillator that initially lies in a thermal state 
with temperature f3o- Then at t ~ we quench its fre- 
quency from uo to uj. The time evolution of the posi- 
tion operator after the quench is given by the Heisenberg 
equations of motion x{t) +uj'^x(t) = with solution x{t) — 
x{0) cos ujt+p{0)smujt / Lo . The computation of the time or- 
dered correlator T{x{ti)x{t2)} thus reduces to finding the 
values of x^, xp, px andp^ in the initial state. These values 
are thermal statistical averages of quantum expectation 
values given in general by 0{f3) = Tr{e-'^"d)/Z{l3) = 



[n\0\n) / Z {(3) . For a simple harmonic oscilla- 



tor of frequency tJo and at temperature /3o, since the en- 
ergy levels are En = {n -\- the partition function is 
Z{f3o) = l/(2sinh/3oWo/2). From the equipartition the- 
orem and the equation E{P) — —d{\nZ{P))/df3 for 
energy of a thermal state, or by direct calculation, 
find 
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Also = i and ({a;,p})/3(, = 0. Using all of the above 
equations we can show that the time ordered correlator is 
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The first term is the Feynman propagator in O+I dimen- 
sions and is the only one that survives if cjq ~ and 
/3o — ^ cxD. The breaking of time-translational invariance 
due to the quench is apparent in the last term that de- 
pends on t\ + 12. Also notice that for /Jq ^ c» we recover 
the pure quantum quench propagator, since the initial 
state is then the ground state of the initial hamiltonian 
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4cc>oC^2 — v^. ' (2) 

Obviously the harmonic oscillator does not thermalize 
but oscillates. However, in a free field theory, we shall 
argue that the interference between momentum modes 
needed to calculate local observables is responsible for 
effective thermalization. 

From the above results we find that the energy of the 
oscillator before the quench is 

(i/o>A, = I^oCoth(/3ot^o/2) (3) 
while after the quench it is 



coth(/3oc^o/2) 



(4) 



i.e. it changes by a factor {E — Eo)/Eo ~ {{uj/ujo)'^ — l)/2. 
This means that the work done on the system at the 
quench is positive or negative, depending on whether the 
frequency increases or decreases respectively. Notice that 
for <C ojo the system loses half of its energy. 

Bosonic free field theory. Let us now move on to 
a bosonic free field theory with hamiltonian of the general 
form 
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and assume a relativistic dispersion relation oJf, ^ (?k^ + 
m?c^ in which we quench the mass from mo to m. We 
will later discuss other types of quenches with different 
dispersion relations or different quench parameters, but 
we now specialize to this case which contains most of the 
physical features of the general cases. For brevity we set 
the speed c = 1. 

The propagator in the mixed momentum-time represen- 
tation is just that of a single harmonic oscillator (1) with 
frequencies ujk and cjofc defined by m and mo respectively. 
The real space propagator is the Fourier transform of the 
latter and we can easily check, using the stationary phase 
method, that for large times and finite separations the in- 
tegration over all momenta leads to the (ii-|-i2)-dependent 
part of the propagator decreasing with time under quite 
general conditions (for m or even for m = in ^). 



^Even though this condition docs not hold for m = 0,rf = 1, we 
can still talk about thermalization in la! critical systems, since what 
is physically meaningful is not the propagator, but the correlation 
function of vertex operators, i.e. imaginary exponentials of the field 
operator, and the latter does become stationary [2]. 



Quantum Quench from a Thermal Initial State. 



This property is already known for the pure quench prop- 
agator (2) [1,2, 14] and we now simply observe that it also 
holds in the more general thermal case (1). According to 
the above, the correlation function tends for large times 
to the stationary form 
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Furthermore one can compare this to the thermal or Mat- 
subara propagator in real time which gives the correlations 
in a system at thermal equilibrium 
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Again as in the pure quench case, the stationary prop- 
agator (6) is of thermal form, where the effective tem- 
perature PcS is determined by equating the coefficients of 
cosujk{ti — t2) in the compared expressions, giving 

r. {i0k-^ok)^+e''"^o-{u;k+ujok)^ 

/^eff ifej - m ^ ^^^^^ ^ ^^^^^^ _ ^^^^^ . [6) 

The effective temperature obtained in this way turns out 
to be a function of the momentum: in a free theory the dif- 
ferent momentum modes do not interact with each other 
and there is therefore no reason why they should all ac- 
quire the same temperature. As before this property is not 
new; it is also true for a pure quantum quench where [2] 
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Note that the three temperatures /3o7/3off and f3*g satisfy 
the following symmetric relation 

/3cff(fc)wfc 



tanh^^i^tanh.^°""'= 



tanh ■ 



(10) 
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Let us investigate the main features of the effective tem- 
perature. We can distinguish the following limiting cases 

• Cold quench wpfe ^ Po^' asymptotically repro- 
duce the pure quantum quench result (9). 



• Hot quench /3q ^ cjok'- 



2/3o 



(11) 



l + ^fcKfe' 

For ujQk > ^k, this leads to the interesting conclusion 
that the system becomes 'colder' after the quench, un- 
like the pure quench case where the system becomes 
'hotter' since the temperature rises from zero to a 
finite value. In particular in the deep quench limit 
ijjQk 3> Wk the final temperature is half of the initial 
one. This can be explained by the comments after 
(4) along with the fact that from (3) the energy of a 
harmonic oscillator at thermal equilibrium with tem- 
perature /3~^ ^ oj tends to the classical value f3~^. 




Figure 1: Effective temperature after a thermal quantum 
quench as a function of took for several values of Po at lo^ = 1. 
Inset: Pcs/Po always for tOk ~ 1. 



• Critical evolution to = 0: In this case, the zero- 
momentum effective temperature has the simple form 

/?cff(fc = 0) = — tanh^, (12) 
Too 2 

showing explicitly the crossover between cold and 
hot quench. In particular, when (S^rriQ = 3.83002 . . . 
the effective temperature is equal to the initial one. 

We point out that the relation jS^s ~ 2/3o for the hot 
critical quench, is the same obtained (for any observable) 
from the classical fluctuation-dissipation relation in the 
ageing regime in a gaussian field theory [17]. Fig. 1 shows 
plots of /Soff as a function of the frequency wofc for sev- 
eral values of /3o. Note that the assumption wofe ^ ^k 
requires both m rriQ and k <^ toq. Although the small 
momentum modes indeed play the most important role in 
determining the stationary behaviour of the system, one 
has to check the contribution of large momenta as well 
and this is taken into account in the following calculation. 

An alternative effective temperature. Ideally 
the effective temperature should be a single number de- 
scribing the asymptotic state after the quench. It is then 
important that this number should not depend on the par- 
ticular quantity we measure. The zero- momentum /3eff(fc) 
is by definition in a free theory a good measure for the 
large-distance behaviour of any correlation function (we 
have derived it from the two-point function, but in a free 
theory this is enough to calculate all of them). Also in a 
conformal field theory the zero-momentum effective tem- 
perature is well-defined and independently observable [2]. 
However, one could wonder whether physical quantities 
involving higher momentum modes can spoil this result. 
To this end, we propose another estimate of the effective 
temperature by the comparison of the field fluctuations 
(0^(x = 0, t ^ oo)) in the two cases: after the quench and 
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Table 1; Asymptotic behaviour of the integrals /i, I2 and Ith ^n 
the limit (3^^ ^ mo S> m. 
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uokieP'^'- - 1) 

Sueh an estimate involves an average over all momentum 
modes. For this reason we will call it 'average' efFcetive 
temperature and denote it by (3. In the massless limit, 
both integrals are infrared divergent in Id or 2d, but the 
two sides compensate each other. 

Eq. (13) can be solved numerically or even analytically 
in the interesting limit /J,^^ ^ rriQ ^ m mentioned above. 
Then we expect f3 to be small too and we can use the 
relevant asymptotic form of the thermal integral 

J _ r k'^-^dk 
"'^Jo c.,(e'5-'^-l)' 

which is easily worked out for any number of spatial di- 
mensions d (see table 1). The integral of the left hand side 
of (13) can be split into two parts: the pure quench part 
which can be calculated exactly 
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and the remaining /3o-dependent part which turns out to 
dominate in the above limit 



k'^-^dk 
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coth^-1 



Table 1 summarizes the asymptotic behaviour of the above 
integrals in the hot deep quench limit. The comparison 
shows that in Id and 2d, (3 = 2/3o as happens for /?cff(fc) 
for small momenta. In 3(i however, (3 = (3q i.e. the average 
temperature tends to the initial temperature instead of its 
half. As a conclusion we observe that, even though the ef- 
fective temperature for small momenta is always half of 
the initial one, the averaging of momenta involved in the 
calculation of the field fluctuations and the fact that in 3c? 
the dominant role of small momenta is reduced, result in 
the average effective temperature remaining the same as 
the initial one. Note also that in 2d there are logarith- 
mic corrections to the asymptotic behaviour which could 
render comparison with data difficult. 

Interestingly enough, in the above limit (3 = /?off(fc = 0) 
in Id and 2d (and, as we will see, in any dimension for a 
lattice model for a sufficiently deep quench). This suggests 
that (3eS could have a sounding physical meaning even 
beyond the small momentum features. 



Two-point correlation in real space. It is worth 
to give a quick look at the equal-time correlation in real 
space. The reason is twofold. Firstly we can explicitly see 
how the zero- momentum effective temperature Pcsik = 0) 
describes its large distance asymptotic. Second, we can 
understand how the pure quench horizon effect [1,2,4] is 
smoothed by finite temperature. 

From (1) the real space propagator for ti = ^2 = ^ is 



d'^k 



ik-r'^k 



+ + (to^ — TOq) cos 2(jjfct 



coth 



/3o^ofe 



i2nY Aivliuok 

where, compared to the pure quench [2] , there is only the 
additional coth factor. Let us firstly consider the massless 
(conformal) evolution with m = and d = 1. For large 
r,t and t > r/2 a saddle point argument gives that the 
correlation function is linear with a slope given by /3eff(0) 
as C{r,t) ^ /3off(0)~ {t — r/2). Fig. 2 shows how this limit 
is reached for several quenches. Notice that in this regime 
the full correlation function is described only by /3off(0), 
that is a single number encoding all features of the quench. 
The massive evolution gives a less trivial saddle point that 
superimposes slowly decaying large oscillations that are 
analogous to those found in the pure quench [2] and will 
not be discussed further. 
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Figure 2: Smoothing of the horizon for real-space correlation. 
All curves refer to critical evolution m = and d = 1. The 
straight lines are the asymptotic behaviour. 

From Fig. 2 it is also evident how the finite temperature 
smooths the horizon, which is not sharp as in the pure 
case. The interpretation of this fact is straightforward: 
quasi-particles emitted at i = from a distance smaller 
than the thermal correlation length cx Pq are entangled 
and generate correlations between two points at distance 
r faster than particles emitted from the same point (or 
inside the true correlation length ttiq^ <^ f3a). One is then 
tempted to conclude that the effect of the temperature is 
similar to ttiq^, but this is only partially true. In fact a 
finite mo, not only smooths the horizon, but also produces 
a shift in the zero time correlation (see the corresponding 
curve in the figure). 

Free bosons on the lattice. — Let us now consider 
a system of free bosons on a lattice whose hamiltonian, if 
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expressed in terms of the field operators, takes the form 

h=j: (ini+i^ui). (14) 

keBZ ^ ^ 
with a lattice dispersion relation, typically of the form 

^2 

Lul = m^c^ + 2—'^{l-coskja), (15) 

j 

where m is again the energy gap, a is the lattice spacing 
and j enumerates the d space coordinates. For convenience 
we set c = 1 and a = 1. 

The only essential difference with the previous treat- 
ment is in the dispersion relation and the fact that the 
momentum k now takes values in the first Brillouin zone 
(BZ) [— TT, +7r]'^. This obviously affects all quantities that 
involve a sum over all momenta, but not the small mo- 
mentum ones. First of all we have to verify that the 
(ii +i2)-dependent part of the propagator is still decaying 
with time. Although there are now additional stationary 
points in cuk that determine the asymptotic behaviour of 
the fc-integral (specifically the edges of the Brillouin zone) 
their contributions also decay with time as before. There- 
fore thermalization also occurs in the lattice. Secondly the 
computation of the average effective temperature is easier 
since the natural cutoff introduced by the lattice spacing 
simplifies the asymptotic analysis of the momentum inte- 
grals. More explicitly, in the deep quench limit, we can 
assume that m-o is much larger than all k in the Brillouin 
zone. Then the left hand side of (13) can be written in 
any dimension as (mo/4) coth(/3oTOo/2) X]fc(l/'^fe) ^^'^ 
right hand side becomes in the same limit J2ki'^/^k)- 
Thus the average effective temperature is 

= imoCoth(/3omo/2), (16) 

which is the same as (12), showing that in the lattice the 
two effective temperatures are equal in any dimension. 

Other quenches. — Up to now we have considered 
only quenches of the mass parameter of a relativistic 
dispersion relation. Other possibilities can also be studied 
in exactly the same way using (1) as the starting point, 
checking the thermalization condition and estimating 
the effective temperature in various limits. For example 
one could quench the speed of sound c or assume a clas- 
sical dispersion relation Wfc = A -I- k^ /2m* and quench 
independently the energy gap A or the effective mass 
m*. However not all of these possibilities are physically 
meaningful in the context of a continuous field theory 
where the small scale degrees of freedom have been com- 
pletely ignored as presumably not playing a crucial role. 
This means that the quench should not alter the large-fc 
behaviour of the spectrum. This explains the emergence 
of ultraviolet divergences in the field fluctuations, in any 
dimension, when one considers quenches of the speed 
of sound c. Similar divergences occur at a quench of 



the effective mass m* in the classical dispersion relation 
above, except for Id and A ^ where we find that the 
system thermalizes. This also happens for all dimensions 
in the case of a quench of the energy gap A only. 

The conclusion is that as long as the quench docs not 
alter the large-fc behaviour of the spectrum in such a way 
as to cause divergences, thermalization occurs for a broad 
class of gapped and gapless dispersion relations. However 
the two simple types of dispersion relations we considered 
so far are not sufficiently adjustable to allow for any 
quench that would not affect the large-A: behaviour, other 
than a quench of the energy gap. Conversely, in lattice 
models, because of the presence of a natural momentum 
cutoff preventing divergences, more elaborate quenches 
are possible and well-defined. If for example we perform 
a (pure) quench of the speed of sound of the lattice 
dispersion relation uj^ = A"^ + 2(? ~ cosfcj), with 

A kept fixed, then the propagator thermalizes as 
before, but the momentum distribution of quasiparticles 
as given by the stationary part of the propagator exhibits 
a maximum at a nonzero k and it vanishes at fc = 0. This 
means that after such a quench the initial pure state feeds 
higher excitation levels rather than the lowest ones which 
remain empty. The case A = 0, d = 1 is again special: the 
propagator does not tend to a constant value, i.e. it does 
not thermalize in the usual sense, but rather increases 
for large times as ln(i/a). However, as mentioned also 
earlier, in the gapless regime the correlation functions of 
vertex operators are physically more meaningful^. 

More complex quenches are possible in free lattice mod- 
els with non-trivial spectra. In any case the procedure for 
the computation of the time evolution is the same and gen- 
eral expressions can be readily obtained. Let us consider a 
system of free bosons described by a general hamiltonian 
of the form 

^ = H + \Bu{a^_^a\ -\- a^kttk), (17) 

k 

where Ak , Bk are even functions of k so that H is her- 
mitian. Suppose that the initial temperature is /3o when 
we quench the parameters of the model from > to 
A, B and that we wish to find the time evolution of the 
momentum distribution nk{t) = {a^ak). To this end 
we need to diagonalisc each of the two hamiltonians by 
means of a Bogoliubov transformation. Then to calculate 
Hkit) = Ti{e-'^»"°e+'"*alake-'"*}/Z{Po) we first have 
to express aj^Uk in the basis in which H is diagonal so 
that the action of the time evolution operator exp(— i_fft) 
can be worked out easily. Next we have to rewrite the 
operators in the basis in which Hq is diagonal so that the 
expectation value on the initial thermal state can also be 
worked out, using the fact that the momentum distribu- 
tion in this basis is equal to the Bose-Einstcin distribution 

^In fact, in the Luttingor model [9], the logarithmic behaviour 
of the boson propagator leads to a power-law time decay of the 
fermionic two-point correlation function. 
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coth 



1). We finally obtain 



^ sin^ Ekt 
El 



(18) 



As usual, the momentum distribution does not equilibrate 
itself, but since any measurement is done on a finite region 
of real space, the oscillating part of Uk (t) is unobservable 
for sufficiently large times due to the momentum averag- 
ing that comes into play. This argument however relies 
on the functional form of the excitation spectrum Ek and 
especially its stationary points. If this is such that the os- 
cillating part can be neglected, the observable distribution 
is the stationary part 



nk 



2EokEl 



coth 



Ok 



(19) 



If we had a system of fermions instead of bosons then we 
could use the same form (17) for the hamiltonian but with 
Bk being an odd function of k and imaginary. For clarity 
let us redefine it by the substitution Bk — > iBk so that it 
is real. The energy spectrum would then be El - 
and following the same procedure we would find 



Bl 



'^k(t) = t;— tanh — - — 
^ofe ^ 



—I sin^ Ekt - - 
El 2 



and 



nk = 



2EokEl 



tanh 



Ok 



t;, (20) 



(21) 



The above expressions can be used to determine the mo- 
mentum distribution after any quench of a free lattice 
model. Note that these results are in agreement with ear- 
lier studies (for example [8, 10]). 

Conclusions. — We have explored various cases of 
quenches in free models that lead to stationary behaviour 
for large times and investigated its effective thermal prop- 
erties, also suggesting a way to estimate an average effec- 
tive temperature from the field fluctuations. We showed 
that if the initial temperature of the system is nonzero and 
higher than other parameters then the quench can lower 
it significantly. We expect that several of our findings are 
not only valid for systems admitting a free particle rep- 
resentation, but are general aspects of thermal quantum 
quenches. 
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